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Abstract 

The aim of this paper is to study the d-dimensional stochastic heat 
equation with a multiplicative Gaussian noise which is white in space 
and it has the covariance of a fractional Brownian motion with Hurst 
parameter H G (0, 1) in time. Two types of equations are considered. 
First we consider the equation in the Ito-Skorohod sense, and later in 
the Stratonovich sense. An explicit chaos development for the solution is 
obtained. On the other hand, the moments of the solution are expressed 
in terms of the exponential moments of some weighted intersection local 
time of the Brownian motion. 



1 Introduction 

This paper deals with the d-dimensional stochastic heat equation 

du 1 . d 2 W H ,„ . 

driven by a Gaussian noise W H which is a white noise in the spatial variable 
and a fractional Brownian motion with Hurst parameter H G (0, 1) in the time 
variable (see (|2.1|) in the next section for a precise definition of this noise). The 
initial condition uq is a bounded continuous function on R d , and the solution will 
be a random field {ut. x , t > 0, x G M. d }. The symbol o in Equation Ijl.ip denotes 
the Wick product. For H = i, 9 d ^ )x is a space-time white noise, and in this 
case, Equation (jl.ip coincides with the stochastic heat equation considered by 
Walsh (see [IT). We know that in this case the solution exists only in dimension 
one (d = 1). 

There has been some recent interest in studying stochastic partial differential 
equations driven by a fractional noise. Linear stochastic evolution equations in 
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a Hilbert space driven by an additive cylindrical fBm with Hurst parameter H 
were studied by Duncan et al. in [3] in the case H £ (|, 1) and by Tindel et al. 
in [T5] in the general case, where they provide necessary and sufficient conditions 
for the existence and uniqueness of an evolution solution. In particular, the 
heat equation 

du_l A d 2 W H 
~dt~2 U+ dtdx 

on K d has a unique solution if and only if H > 4 . The same result holds when 
one adds to the above equation a nonlinearity of the form b(t,x,u), where b 
satisfies the usual linear growth and Lipschitz conditions in the variable it, uni- 
formly with respect to (t, x) (see Maslowski and Nualart in [9]). The stochastic 
heat equation on [0, oo) x ~R d with a multiplicative fractional white noise of 
Hurst parameter H = (Ho, Hi, ... , Hj) has been studied by Hu in [6] under the 
conditions j- < Hi < 1 for i = 0, . . . , d and J2i=o Hi < d — 2 h -i ' 

The main purpose of this paper is to find conditions on H and d for the 
solution to Equation to exist as a real-valued stochastic process, and to 

relate the moments of the solution to the exponential moments of weighted in- 
tersection local times. This relation is based on Feynman-Kac's formula applied 
to a regularization of Equation In order to illustrate this fact, consider 

the particular case d = 1 and H = h. It is known that there is no Feynman- 
Kac's formula for the solution of the one-dimensional stochastic heat equation 
driven by a space-time white noise. Nevertheless, using an approximation of the 
solution by regularizing the noise we can establish the following formula for the 
moments: 



k I k t 

f[u (x + Bi)e^i £ f S (Bl-Bi 

3 = 1 \i,j=l,i<j J ° 



)ds 



, (1-2) 



for all k > 2, where B t is a fc-dimensional Brownian motion independent of the 
spaced-time white noise W?. In the case H > i and d > 1, a similar formula 
holds but J* 5q(BI — B 3 s )ds has to be replaced by the weighted intersection 
local time 



L t = H(2H - 1) [ [ \s - r\ 2H ~ 2 5 Q {B\ ~ B^.)dsdt, 
Jo Jo 



(1.3) 



where {B^ ,j > l} are independent d-dimensional Brownian motions (see The- 
orem 15. 3p . 

The solution of Equation (|1.1[) has a formal Wiener chaos expansion ut, x — 
Y^=o^ n (f n (''^^ x ))- Then, for the existence of a real-valued square integrable 
solution we need 

oo 

^n! \\fn{',t,x)\\ n 8n < oo, (1.4) 

n=0 

where Hd is the Hilbert space associated with the covariance of the noise W . It 
turns out that, if H > i, the asymptotic behavior of the norms ||/ n ( - , t, x)\\ H ®n 
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is similar to the behavior of the nth moment of the random variable L t defined 
in (|1.3p . More precisely, if uq is a constant K, for all n > 1 we have 

(n\) 2 \\f n (-,t : x)\\ 2 nT =K 2 E(U?). 

These facts leads to the following results: 

i) If d = 1 and H > i, the series (|1.4I) converges, and there exists a solution 
to Equation which has moments of all orders that can be expressed in 
terms of the exponential moments of the weighted intersection local times 
L t . In the case H = i we just need the local time of a one-dimensional 
standard Brownian motion (see (|1.2p ). 

ii) If H > i and d < 4H, the norms \\f n (-,t,x)\\ n ®n. are finite and E(L™) < 
oo for all n. In the particular case d = 2, the series (ll.4[) converges if t is 
small enough, and the solution exists in a small time interval. Similarly, 
if d = 2 the random variable L t satisfies -E(expAL t ) < oo if A and t are 
small enough. 

iii) If d = 1 and | < H < ^, the norms ||/n(-, t, x)\\ n ®n are finite and E(L™) < 
co for all n. 



A natural problem is to investigate what happens if we replace the Wick 
product by the ordinary product in Equation that is, we consider the 

equation 

du 1 . d 2 W H 

In terms of the mild formulation, the Wick product leads to the use of Ito- 
Skorohod stochastic integrals, whereas the ordinary product requires the use 
of Stratonovich integrals. For this reason, if we use the ordinary product we 
must assume d = 1 and H > h. In this case we show that the solution exists 
and its moments can be computed in terms of exponential moments of weighted 
intersection local times and weighted self-intersection local times in the case 

The paper is organized as follows. Section 2 contains some preliminaries on 
the fractional noise W H and the Skorohod integral with respect to it. In Section 
3 we present the results on the moments of the weighted intersection local times 
assuming H > i. Section 4 is devoted to study the Wiener chaos expansion 
of the solution to Equation The case if < ^ is more involved because 

it requires the use of fractional derivatives. We show here that if | < H < |, 
the norms t, x)\\~,®n are finite and they are related to the moments of a 

fractional derivative of the intersection local time. We derive the formulas for 
the moments of the solution in the case H > ^ in Section 5. Finally, Section 6 
deals with equations defined using ordinary product and Stratonovich integrals. 
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2 Preliminaries 

Suppose that W H = {W H (t,A),t >0,ie B{R d ),\A\ < oo}, where B(R d ) is 
the Borel cr-algebra of R d , is a zero mean Gaussian family of random variables 
with the covariance function 

E(W H {t, A)W H (s, B)) = ^(t 2H + s 2H -\t- s\ 2H )\A n B\, (2.1) 

defined in a complete probability space (CI, J 7 , P), where H £ (0,1), and \A\ 
denotes the Lebesgue measure of A. Thus, for each Borel set A with finite 
Lebesgue measure, {W H (t,A),t > 0} is a fractional Brownian motion (fBm) 
with Hurst parameter H and variance t \A\, and the fractional Brownian 
motions corresponding to disjoint sets are independent. 

Then, the multiplicative noise d d ^ x appearing in Equation (jl.ip is the for- 
mal derivative of the random measure W (t, A): 

rr r /■' d 2 w H 

We know that there is an integral representation of the form 
W H (t,A)= f f K H (t,s)W(ds,dx), 

JO J A 

where W is a space-time white noise, and the square integrable kernel Kh is 
given by 

K H (t,s) = c H s^~ H (u- s) H -iu H -^du, 

J s 

for some constant ch ■ We will set ifjy(i, s) = if s > t. 

Denote by £ the space of step functions on R + . Let Tt be the closure of £ 
with respect to the inner product induced by 

(![o,t], 1 [o, s ]) n = K H (t,s). 

The operator K* H : £ — * L 2 (R + ) defined by K H (l[ ^)(s) = K H (t,s) provides 
a linear isometry between Ti and L 2 (K + ). 

The mapping l[o.t]xA ~~ * W (t,A) extends to a linear isometry between the 
tensor product H (8) i 2 (R d ), denoted by Hd, and the Gaussian space spanned 
by W H . We will denote this isometry by W H . Then, for each <p 6 Hd we have 

W H (tp)= f f (K* H <E>I)ip(t,x)W(dt,dx). 
Jo Jm d 

We will make use of the notation W H (if) = J °° J Rd ipdW H . 

If H = i, then Ti = i 2 (K + ), and the operator K* H is the identity. In this 
case, we have H d = L 2 (R + x R d ). 
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Suppose now that H > i. The operator K* H can be expressed as a fractional 
integral operator composed with power functions (see [llj). More precisely, for 
any function <p £ £ with support included in the time interval [0, T] we have 

(K* H <p) (*) = c' H ti- H I^ (t), 

H— - 

where I T _ 2 is the right-sided fractional integral operator defined by 

J T- _I /W = j\* ~ t) H ^f(s)ds. 

In this case the space H. is not a space of functions (see [2]) because it contains 
distributions. Denote by \H\ the space of measurable functions on [0, T] such 
that 

oo />oo 

/ \r — u\ 2H ~ 2 \(p r \\ip u \drdu < oo. 
10 Jo 

Then, \H\ C Ti. and the inner product in the space TL can be expressed in the 
following form for tp,ip £ \7i\ 

I* OO />OC 

{ ( P>' l P)n = / / (f>(r.,u)(p r (p u drdu, (2.2) 
Jo Jo 

where <£(s,i) = H {2H - l)\t - s\ 2H ~ 2 . 

Using Holder and Hardy-Littlewood inequalities, one can show (see [TO] ) 
that 

IMIw d </MMU (E+;i2(E<i)) , (2-3) 
and this easily implies that 

ii^ll«r^^lkli^ (R „ iW) - (2.4) 

If H < i, the operator K* H can be expressed as a fractional derivative 
operator composed with power functions (see [11]). More precisely, for any 
function ip £ £ with support included in the time interval [0, T] we have 

[K* H <p) (t) = c'^- H DY_ H (v(s)s H -i) (t), 
where is the right-sided fractional derivative operator defined by 

(r-t)i- H l 2 ' J t ( s - t ) H -i 



Moreover, for any 7 > | — H and any T > we have C 7 ([0,T]) C H 

4: h (l 2 ([o,t}). 
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If if is a function with support on [0, T], we can express the operator Kjj in 
the following form 

K*Mt)=K H (T,t)<p(t) + J [<p( S )-v(t)]^(s,t)ds. (2.5) 
We are going to use the following notation for the operator K* H : 

K* H p= f <p(t)K* H (dt,r). (2.6) 

J [O.T] 

Notice that if H > i, the kernel Kh vanishes at the diagonal and we have 
K* H (dt,r) = 9|»(i,r)l [r>T] (i)dt. 

Let us now present some preliminaries on the Skorohod integral and the Wick 
product. The nth Wiener chaos, denoted by H„, is defined as the closed linear 
span of the random variables of the form H n (W H ((f)), where ip is an element of 
Tld with norm one and H n is the nth Hermite polynomial. We denote by I n the 
linear isometry between H.® n (equipped with the modified norm \fn\ ||-|| w ®n) 

and the nth Wiener chaos H„, given by I n ((p® n ) — n\H n (W H (<p)), for any 
ip G Ti-d with ||y||^ = 1. Any square integrable random variable, which is 
measurable with respect to the a- field generated by W H , has an orthogonal 
Wiener chaos expansion of the form 

oo 

F = E(F) + J2Wn), 

n=l 

where /„ are symmetric elements of TC® n , uniquely determined by F. 

Consider a random field u = {u t<x ,t > 0, x £ R d } such that E (m| x ) < oo 
for all t, x. Then, u has a Wiener chaos expansion of the form 

oo 

u t , x = E{u t . x ) + In(fn(; t,x)), (2.7) 

n=l 

where the series converges in L 2 (Vl). 

Definition 2.1 We say the random field u satisfying is Skorohod inte- 

grable if E(u) € Tid, for all n > 1, /„ G %^ n+1 \ and the series 

oo 

W H (E(U)) + J2 In + X(fn) 
n=l 

converges in L 2 (Q), where f n denotes the symmetrization of f n .We will denote 
the sum of this series by 8(u) = J Q J Rd u5W H . 
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The Skorohod integral coincides with the adjoint of the derivative operator. 
That is, if we define the space O 1,2 as the closure of the set of smooth and 
cylindrical random variables of the form 

F = f(W H (h 1 ),...,W H (h n )), 

hi G Hd, f G C^°(IR n ) (/ and all its partial derivatives have polynomial growth) 
under the norm 

\\DF\\ 1<2 = JE(F*)+E(\\DF\\ 2 nd ), 



where 

n 

dx, 
3=1 J 

then, the following duality formula holds 



df = J2 ^(W H (hi), . . . , W B {h n ))h j , 



E(6(u)F)=E((DF,u) Hd ), (2.8) 

for any F G D 1 ' 2 and any Skorohod integrable process u. 

If F G D 1,2 and h is a function which belongs to Hd, then Fh is Skorohod 
integrable and, by definition, the Wick product equals to the Skorohod integral 
of Fh: 

S(Fh) = FoW H {h). (2.9) 
This formula justifies the use of the Wick product in the formulation of Equation 

Finally, let us remark that in the case H = i, if Ut, x is an adapted stochastic 
process such that E (J °° J Rd u 2 x dxdt) < oo, then u is Skorohod integrable and 
6(u) coincides with the Ito stochastic integral: 



5{u)= / / u t , x W{dt,dx). 



3 Weighted intersection local times for standard 
Brownian motions 

In this section we will introduce different kinds of weighted intersection local 
times which are relevant in computing the moments of the solutions of stochastic 
heat equations with multiplicative fractional noise. 

Suppose first that B 1 and B 2 are independent <i-dimensional standard Brow- 
nian motions. Consider a nonnegative measurable function f](s,t) on R 2 ^. We 
are interested in the weighted intersection local time formally defined by 



1= [ [ r](s,t)S (Bl - B 2 )dsdt. (3.1) 
Jo Jo 

We will make use of the following conditions on the weight rj: 
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CI) For all T > 



1 T := max sup 

\ o<t<r Jo 



rj(s,t)ds, sup / r)(s,t)dt\ < oo. 

0<s<TJo J 



C2) For all T > there exist constants 7t > and 6 (0, 1) such that 

V(s,t) < It \s - t\ 2H ~ 2 , 

for all s,t<T. 

Clearly, C2) is stronger than CI). We will denote by pt(x) the d-dimensional 

d w 2 

heat kernel pt{x) = (2nt)~ 2 e~^*~ . Consider the approximation of the inter- 
section local time (|3.1[) defined by 

i-T t-T 



I s = [ f v (s,t)p e (Bl - B 2 )dsdt. 
Jo Jo 



(3.2) 



Let us compute the fcth moment of I e , where k > 1 is an integer. We can write 

k 



E ( I e)= I f[r)(si,ti)ilJ e (s,t)dsdt, 

J]0.T] 2k 



where s = (si, . . . , s k ), t = (ti, . . . , t k ) and 

lk(s,t) = E {p^B^ - Bl) ■ ■ - Pe (Bl k - B 2 J) 
Using the Fourier transform of the heat kernel we can write 



(3.3) 



(3.4) 



V>e(s,t) 



(27T) 



kd 



R k 



e ex p £K^~^)~f ^ 



(2 



7r) \JR k 



(3.5) 



where £ = (£1, • • • , £fc) an d b't, i = 1,2, are independent one-dimensional Brow- 
nian motions. Then ip £ (s,t) < V>(s,t), where 



■0(s,t) = (2tt) _ ~ [det (sj A si+tj A U)} 2 



Set 



Oik 



\\r){s i ,t i )ip{s,t)dsdt. 

[o,T] 2 * l=1 



(3.6) 



(3.7) 



Then, if a k < 00 for all k > 1, the family J e converges in L p , for all p > 2, to a 
limit / and E(I k ) — a k - In fact, 

lim E(I e Ig) = a 2 , 

so I e converges in L 2 , and the convergence in L p follows from the boundedness 
in L q for q > p. Then the following result holds. 
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Proposition 3.1 Suppose that CI) holds and d = 1. Then, for all A > the 

random variable defined in \3. || ) satisfies 



supB(exp(AJ e ))<l + $^|H| 1)T Aj , (3.8) 

where Q(x) — ^2kLi r ^+i ^ ■ -<4iso, 7 e converges in LP for all p > 2, and the 
limit, denoted by I, satisfies the estimate Ii3.8\) . 

Proof The term ip(s,t) defined in (|3.6p can be estimated using Cauchy- 
Schwarz inequality: 

ip(s,t) < (2ir)~ k [det (sj A s;)]^ [det (tj A ij)]* 

= 2- fe 7r-i[/3(s)/3(t)]-i, (3.9) 

where for any element (si, . . . , s&) € (0,oo) fe with s, 7^ Sj if i 7^ j, we denote 
by cr the permutation of its coordinates such that s a m < ■ ■ ■ < s a ^ and 
(3(s) = s (T (i)(s (T (2) - s^x)) • • ■ (s cr(fc ) - s (T(fc _ 1) ). Therefore, from (|3~9l) and (|3~7| 
we obtain 

a k <2- k n-i [ TT 77(5,, fc) [/3(s)/?(t)r* dsdt. (3.10) 
Applying again Cauchy-Schwarz inequality yields 

a fe < 2- k n-U [ T\r)( Si ,ti)[l3(s)]-i dsdt] 

TT^s^H/^Wr^sdtl 

[0,Tp J 



< (2-^-^ \\n\\ 1T ) k kl f [/3(s)]-»dB 

rtfj ' ' 



(3.11) 



where Tj. = {s = (s\, . .., s&) ■ < Si < • • • < < T}, which implies the 
estimate (|3.8[) . ■ 

This result can be extended to the case of a (i-dimensional Brownian motion 
under the stronger condition C2): 

Proposition 3.2 Suppose that C2) holds and d < AH . Then, lim £ |o le = I, 
exists in LP , for all p > 2. Moreover, if d — 2 and A < \q(T), where 

, , , H(2H-l)An d 9W / d V 2H , _ , 

A "^= 7r/3g r? F ( iff) ^ (3 ' 12) 
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and (3h is the constant appearing in the inequality i2.3\) , then 

sup£(exp(A/ £ )) < oo, (3.13) 

and I satisfies E (exp (XI)) < oo. 

Proof As in the proof of Proposition 13.11 using condition C2) and inequality 
(|2.4p we obtain the estimates 

a k < 7 £2- dfe 7r^ / T\\t t - Sl \ 2H - 2 [f3( S )P(t)}-Usdt 

J[0,T]« ZJ[ 

< ^2- dk n-fa k H ( [ [/3(s)]-^ds] 

\J[0,T]>° J 



( lT a H 2- 2 7r- 1 ) k (fc!) 



k ^ ^ T{1 _^ )k 2H T k(l-^)2H 

T(k(l-^L) + l)^ 



k 



where a^r = h ^h-i) • This allows us to conclude the proof. ■ 

Ifd = 2and?7(s,i) = 1 it is known that the intersection local time J Q T Jq 8q(BI- 
B 2 )dsdt exists and it has finite exponential moments up to a critical exponent 
Ao (see Le Gall [7] and Bass and Chen PQ). 

Consider now a one-dimensional standard Brownian motion B, and the 
weighted self-intersection local time 

i-T 



As before, set 



I=[ [ n(s,t)S {B s ~ B t )dsdt. 

Jo Jo 

I e = f [ V {s,t)p E (B s - B t )dsdt. 

Jo Jo 



Proposition 3.3 Suppose that C2) holds. If H > |? then we have 

sup E (exp (A [J e - E (/,)])) < oo, (3.14) 

£>0 

for all A > 0. Moreover, the normalized local time I — E (I) exists as a limit in 
L p of I £ — E (l e ), for all p > 2, and it has exponential moments of all orders. 
If H > j, then we have for all A > 

sup-E 1 (exp (Alg)) < oo, (3.15) 

for all A > 0, and the local time I exists as a limit in L p of I e , for all p > 2, 
and it is exponentially integrable. 



10 



Proof We will follow the ideas of Le Gall in [7] . Suppose first that H > \ and 
let us show (|3.14p . To simplify the proof we assume T = 1. It suffices to show 
these results for 



J e := [ f r)(s,t)p e (B s - B t )dsdt. 
Jo Jo 



Denote, for n > 1, and 1 < k < 

"2fc-2 2k- 1 



A, 





'2k- 1 2k' 


X 






2 n ' 2™ 



Set 



and 



X s k = I tj{s, t)p e (B s - B t )dsdt 
Ja„ k 



<,& = «n,fc ~ ^ («n,fc) • 

Notice that the random variables a^ fc , 1 < k < 2™~ 1 , are independent. We 
have 

oo 2™ _1 
n=l k=l 

and 

oo 2™ _1 

j e -s(j e ) = X;x;^.fc- 

n=l fc=l 

We can write 



n.k 



= 2 



-2)1 



/ ^ \ 9™ 9™ ' 9^ 

JO \ z z z 

y.y F {B ik-i — Bzk-i, t )dsdt 

on on on "+" 



/2fc- 1 



2fc - 1 i 



2" 



< 7l 2- 2n -( 2H ~ 2 )™ 



which has the same distribution as 



B ik-\ , t )dsdt, 

on +on-' 



1 1 \ t+ ,r~*MB^ 

>n as 

" f ( T \t + s\ 2H - 2 Pt2 n (Bl - B 2 )dsdt, 
Jo Jo 



where B l and B 2 are independent one-dimensional Brownian motions. Hence, 
using the estimate (13.111) . we obtain 



£(exp (A = l + E^fK 



i=2 J 



< 



< 



™ (c T 2-i n ^ 2H - 2 >\Y 

i + E 1 ^ — L 

J=2 
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for some constant Ct- Hence, 

E (exp (A (fi^ fc ))) < 1 + CA 2- 3 "- 2 ( 2ff - 2 )"A 2 , (3.16) 

for some function c\. 

Fix a > such that a < 2{2H - 1). For any N > 2 define 

N 

b N = H(l-2-<i-V), 

and notice that limAr^oo&jv = &oo > 0. Then, by Holder's inequality, for all 
N > 2 we have 



TV 2" 



N-l 2" 



x ^ E 



< IE 



exp A6jv E 

\ 71=1 fe=l 

eX P 1 _ o-Vl) 51 E «n, 
\ n=l fc=l 

exp Xb N 2< N ^ *k 

y fe=i 

/ W-12" -1 

exp A6jv_i ^ a^ fc 



1 _ 2 -a(N-l) 



ra=l fe=l 



{i?[exp (a^^" 1 ^)] } 



2 (l-o)(N-l) 



Using (|3.16p . the second factor in the above expression can be dominated by 

2 (l- a )(N-l) 



{^[exp^&^^W.fc)]} 

< (l + CA A 2 6 2 2 2 < JV - 1 )2- 3JV - 2 ( 2ff - 2 ) JV ) 
(W^- 2 - 2 ^- 2 ^) , 



2 (l-a)(JV-l) 



exp i 



where n — b\2 a 1 . Thus by induction we have 



E 



Xb N Y E < exp \Y KCA A 2 2(- 2 - 2 ( 2ff - 2 ))" 

n=l fe=l y l n=2 J 



x_E (exp ai,x) 
< exp( K c A A 2 (l-2 a+2 - 4ff )- 1 ) 
xE (exp(ai i i)) < oo, 
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because a < 2(2H — 1). By Fatou lemma we see that 



supE (exp (A&oo (J e - E (J £ )))) < oo, 

e>0 



and IpTTiji follows. 

On the other hand, one can easily show that 

Urn E((J £ - E (J £ )) (J 5 - E ( J 4 ))) - i- 

e,oJ.O Z7T 



r)(s,t)r)(s',t' 



s<t<l,s'<t'<l 



det 



t-« |[s,i] n [s',£']| 
|[ s ,t]n[ s ',i']| f-a' 



((i- s )(t'- s '))- 



dsdtds'dt' < oo, 



which implies the convergence of I £ in L 2 . The convergence in L p for p > 2 and 
the estimate Q3.14p follow immediately. 

The proof of the inequality (|3 . 1 5[) is similar. The estimate (|3.16[) is replaced 

by 

E (exp (A « fc ))) < 1 + d A 2- 3 - 2 ( 2ff - 2 )"A, (3.17) 
for a suitable function d\, and we obtain 



E 



N 2" _1 

exp ( Xb N ^2 ^2 an > k 

n=l fe=l 



< exp |^^d A A2(-i- 2ff ) n | E (exp (a M )) 

< exp(V^d A A 2 (l - 2(-^ 2H )")- 1 )£; (exp (a M )) < oo, 
because H > |. By Fatou lemma we see that 

sup -E (exp (Afroo (J e - E (J £ )))) < oo, 

which implies (|3.15[) . The convergence in L p of J £ is proved as usual. 
Notice that condition H > § cannot be improved because 



\t — s\ 1 dsdt = oo. 



4 Stochastic heat equation in the Ito-Skorohod 

sense 

In this section we study the stochastic partial differential equation (|1.1[) on 
R d , where W H is a zero mean Gaussian family of random variables with the 
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covariance function (|2.ip , denned on a complete probability space (n,!F,P), 
and the initial condition uo belongs to Cb(M. d ). First we give the definition of a 
solution using the Skorhohod integral, which corresponds formally to the Wick 
product appearing in Equation (|1.1[) . 

For any t > 0, we denote by Tt the a- field generated by the random variables 
{W(s,A),0 < s <t,Ae B(R d ), \A\ < 00} and the P-null sets. A random field 
u = {ut yX , t > 0, x £ M.} is adapted if for any (t, x), Ut, x is ^"t-measurable. 

For any bounded Borel function tp on K we write pt<f{x) — J Rd pt(x — 

y)<p{y)dy. 

Definition 4.1 An adapted random field u = {ut, x >t > 0,x E K d } such that 
E(iq x ) < 00 for all (t,x) is a solution to Equation hl.l\) if for any (t,x) € 
[0,oo) x R d , the process {pt- s (x — y)u s , y l[o t t](s),s > 0, y S R d } is Skorohod 
integrable, and the following equation holds 

u t ,x — Ptua(x) + / / pt- s {x - y)u StV 6W" y . (4.1) 

JO JM d 

The fact that Equation (jl.ip contains a multiplicative Gaussian noise allows 
us to find recursively an explicit expression for the Wiener chaos expansion of 
the solution. This approach has extensively used in the literature. For instance, 
we refer to the papers by Hu [BJ, Buckdahn and Nualart [2J, Nualart and Zakai 
[13], Nualart and Rozovskii [T2J, and Tudor [TBJ, among others. 

Suppose that u — {u tiX1 t> 0,x & R d } is a solution to Equation (jl.ip . Then, 
for any fixed (t, x), the random variable Ut, x admits the following Wiener chaos 
expansion 

00 

«t,x = £ *«(/«(-, *,*)), (4-2) 

where for each (t, x), /„(•,£, x) is a symmetric element in Hf n . To find the 
explicit form of /„ we substitute (|4.2|) in the Skorohod integral appearing in 
(|4.ip we obtain 

/ Pt-(«-I/)«.,»«W^ = E/V 4(ft- s (i-!/)U,s,!/)X 

oc ^ 

= ^2ln+i{Pt-s(x - y)f n (-,s,y)) . 

71=0 

Here, (p t - s (x — y)f n (', s, y) denotes the symmetrization of the function 
Pt- S (x - y)f n (si,xi; . . . ; s n , x n ; s, y) 
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in the variables (si,xi), . . . , (s n ,x n ), (s,y), that is, 

Pt-s{x - */)/»(•, s, y) = — — -\p t -s(x - y)f n {si,Xi, . . . ,s n , x n , s, y) 
n + 1 

n 
3=1 

x/«(si,xi, . . .jSj-i^j-i^jj/jSj+ijOTj'+i, . . .,s n ,y n ,Sj,yj)]. 
Thus, Equation (|4.1[) is equivalent to say that fo(t,x) = ptUo(x), and 

/ n+ i(-,t,x) = p t -s{x -y)fn{-, s,y) (4.3) 

for all n > 0. Notice that, the adaptability property of the random field u 
implies that f n {si, x\, . . . , s n , x n , t,x) — if Sj > t for some j. 
This leads to the following formula for the kernels /„, for n > 1 

fn ($1: ■ • ■ ) Sn: X n , t , x) - 

rt! 

XPt- M „j (a: - I<r(n)) ' ' f !l(r » l(1] (^tr(2) ~ Z<r(l) «oO CT (l) ) , (4.4) 

where a denotes the permutation of {1, 2, . . . , n} such that < s CT m < • • • < 
s a(n) < £• This implies that there is a unique solution to Equation (|4.ip . and the 
kernels of its chaos expansion are given by (|4.4[) . In order to show the existence 
of a solution, it suffices to check that the kernels defined in (|4.4[) determine 
an adapted random field satisfying the conditions of Definition 14.11 This is 
equivalent to show that for all (t, x) we have 

oo 

^n! \\f n {-,tix)\\ 2 H ®n < oo. (4.5) 

It is easy to show that (|4.5j) holds if if = ^ and rf = 1. In fact, we have, assuming 
|wo| < K, and with the notation x = {x\, . . . ,x n ), and s = (s\, . . . ,s n ): 

\\fn{;t,x)\\ n »n 
1 



{n\f J[ ,t] 



Pt-srw (x ~ x a{n) ) ■ ■ -p SCT(2) _ s<T(1) (x ct(2 ) - ar^i)) 
xp s<T(1) u (a; o -(i)) 2 dxds 



(n!) ./[o,t]»£5 



/T "3 = 1 



where T n — {(si, . . . , s„) € [0, t] n : Q < si < • • • < s n < t} and by convention 
s n+ i = i. Hence, 



\\fn(-,t,x)\\L®n < 



K 2 2~ n t% 



nr - „ir(2±=)' 
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which implies (|4. 5|) . On the other hand, if if = i and d > 2, these norms are 
infinite. 

Notice that if uo = 1, then (n!) 2 ||/ n ( - , t, x)\\^®n coincides with the moment 
of order n of the local time at zero of the one-dimensional Brownian motion 
with variance 2t, that is, 



(n\) 2 \\f n (;t,x)\\ 2 nr =E 



S (B 2s )ds 



To handle the case H > ^ , we need the following technical lemma. 
Lemma 4.2 Set 

g s (xi, x n ) = p t - SaM (x - x a(n) ) ■ ■ •p M2) - Ml) (x a{ 2) - ^(l)))- (4-6) 
Then, 

(Ss,fft) L 2 (R „ d) = V(s,t), 

where ip(s,t) is defined in l{3.4\ )- 
Proof By Plancherel's identity 



(3s,5t) L 2( K „d) = (2vr) ™ (Tg sl Tg t ) L -2 (w , nd) , 
where J 7 denotes the Fourier transform, given by 

n 

J'=l 



n cxp ( * & 

' j=i V 



2 (s CT (j+i) - s ff (j)) 



rfx, 



with the convention = x and s„+i = i. Making the change of variables 



^ ^o-fj) if 1 < J < — 1, and u n — x — x a r n \, we obtain 



•^(Cl, ■ ■ • , £n) = (2tt) "2 JJ(s CT y +1 ) - S ct(j -))" 



/ n ( 



= E I ]^[ exp (t ($ a(j) , a? - B t - )) 



U' =1 



= £ [ J]exp(i(&,a:-B t -.B > ,» 

u=i 
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As a consequence, 

^ftUn^l^l "" / El I I «'X]> ( / U;.Bl : ■ B'f ) ) \ 



which implies the desired result. ■ 

In the case H > ^> and assuming that uq = 1, the next proposition shows 
that the norm (n,!) 2 ||/ n (-, t, x)\\^®n coincides with the nth moment of the in- 
tersection local time of two independent c£-dimensional Brownian motions with 
weight (j>(t, s). 

Proposition 4.3 Suppose that H > A and d < AH. Then, for all n > 1 



(ri) 2 \\f n (;t,x)f H? n < \\uo\LE 



<t>(s,r)6 (Bl - B 2 r )dsdr 



< oo, 
(4.7) 



with equality if uq is constant. Moreover, we have: 

1. If d — 1, there exists a unique solution to Equation 

2. If d — 2 , then there exists a unique solution in an interval [0, T] provided 
T < Tq, where 



d 



-1/(2H-1) 



T =^r(l-— . (4.8) 

Proof We have 

(«0 2 \\fn(-,t,x)\\ 2 H » n < \\uoWlc / YiHsj^j) (g s ,gt) L 2 (Rnil) dsdt, (4.9) 



where g s is defined in (|4.6[) . Then the results follow easily from from Lemma 
[4721 and Proposition [3721 ■ 

In the two-dimensional case and assuming H > A, the solution would exists 
in any interval [0, T] as a distribution in the Watanabe space O"' 2 for any a > 
(see QI]). 

4.1 Case H < \ and o! = 1 

We know that in this case, the norm in the space TL is defined in terms of 
fractional derivatives. The aim of this section is to show that ||/n("> t, x)\\jj»n 
is related to the nth moment of a fractional derivative of the self-intersection 
local time of two independent one-dimensional Brownian motions, and these 
moments are finite for all n > 1, provided I < H < A. 
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Consider the operator (Kg) on functions of two variables denned as the 
action of the operator K H on each coordinate. That is, using the notation (|2.5p 
we have 



(K* H f 2 f(n,r 2 ) = K H (T, n)K H (T, r 2 )/(r 1; r 2 ) 
+K H (T,n) [ ^-(s,r 2 ) (fin, s) - /(n,r 2 ))< 



/"* dK H 

+K H {T,r 2 ) J —^-(v,n)(f(v,r 2 )-f(ri,r 2 ))dv 

1 rt dK H dK H 

— (s,r 2 )— — (v,n) [f(v,s) - f{n,s) - f(v,r 2 ) - f(n,r 2 )\ dsdv. 

r 2 j ri OS OV 

Suppose that f(s, t) is a continuous function on [0, T] 2 . Define the Holder norms 

f \f(s 1: t)-f(s 2 ,t)\ 
II/II17 =sup<^ : ,si,S2,t e T,si ? s 2 

I \ s l ~ s 2\' 
Il/ll2, 7 - 8U P | U _ ,tl,t 2 ,S E l,tl yt t 2 

and 

|/(*i,ti)-/(*i,*2)-/(«a,ti) + /(*2,*2)| 



'1,2,7 



sup ■ 



|*i -s 2 | 7 |ti -i 2 | 7 



where the supremum is taken in the set {ti,t 2 , s 2 , s 2 € T, si ^ s 2 ,t\ t 2 }. Set 

ll/llo, 7 = ll/lll, 7 + ll/lk 7 + ll/lll,2, 7 

Then, (K* H f 2 f is well defined if \\f\\ 0n < oo for some 7 > § - H. As 
a consequence, if -B 1 and B 2 are two independent one-dimensional Brownian 
motions, the following random variable is well defined for all e > 

Je = f {K* H f 2 p £ (B l - B 2 )(r,r)dr. (4.10) 
Jo 

The next theorem asserts that J £ converges in LP for all p > 2 to a fractional 
derivative of the intersection local time of B 1 and B 2 . 

Proposition 4.4 Suppose that | < H < ^.Then, for any integer k > 1 and, 
T > we /iaue £ ( J £ fe ) > and 

sup£ (J e fc ) < 00. 

Moreover, for all p > 2, J 6 converges in LP as e tends to zero to a random 
variable denoted by 

1 (K* H f 2 6o(B} - B 2 )(r,r)dr. 
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Proof Fix k > 1. Let us compute the moment of order k of J e . We can write 
E(j k £ ) = [ E [Y[(K* H f 2 p 2e {B x -B 2 )(r iiri ) ) dr. (4.11) 
Using the expression (|2.6[) for the operator K* H , and the notation (|3.4[) yields 
E(J*)= [ ^ 8 (s,t)TT^(d« 4) r i )A'^(dt i ,r i )tir. (4.12) 

■/[0,T]3* f = i 

As a consequence, using (|3.5[) we obtain 
£ (J £ fc ) - (27r)- fc 



- £*,=! fe«iCov(si . -B t 2 . .B* -B 2 H ] 



J[0,T] k JW k J[0,T] 2k 
k 



< (27T)- 



[0,T]» 



Then, it suffices to show that for each k the following quantity is finite 

/ f 2 \{[s ] -s 3 - 1 +t ] -t^ 1 ]~^WK* H {ds l ,r l )\{K* H {dt l ,r l )dv, (4.13) 

where T k = {0 < t\ < ■ ■ ■ < tk < T}. Fix a constant a > 0. We are going to 
compute 

/ U^-t^+a^YlK^dt^n). 

•' Tk j=l i=i 

To do this we need some notation. Let Aj and 7j be the operators defined on 
a function f(t%, . . . ,t k ) by 

A J -/ = /-/k=r 3 -, 

and 

^/ = I\t j= r r 

The operator K^ 1 {dt il r i ) is the sum of two components (see (|2.5|l ). and it suffices 
to consider only the second one because the first one is easy to control. In this 
way we need to estimate the following term 

/ / Ai-.-Afe [ H tf~* [tj -tj-i +a]-*l {tj _ 1<tj} 

JT k J[0,T]X \ J = 1 



dr. 
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Because tj 2 r? < 1, we can disregard the factors r ? and tj 2 . Using 
the rule 

Aj(FG) = - F{ rj )G{ rj ) 

= [F^) - F( rj )} Git,) + F( rj ) [Git,) - G(rj)] 
= AjFG + IjFAjG, 

we obtain 



s j=i 



where Sj is an operator of the form: 



Ai ■ ■ ■ A fc (jj [tj - tj-x + a] * l {tj _ 1<tj } 
\i=i 

k 

EnM |iri, - i+ " rli fH<' 1 ))' 



II, .IX. A, : I,.X ,A,. 

and for each j, Aj must appear only once in the product Y[j=i Sj- Let us 
estimate each one of the possible four terms. Fix e > such that H — | > 2s. 

1. Term II, r. 

n 3 {\pj - + °]~' i {ti-i<*j}) = fa - + i^j-i^}, 

2. Term /A,-: 



= fe - + a]"' ift-i^} - fa - tj-i + a P 2 1 {t j _ 1 <r j } 
< C - r ]\ 2 ~ H+E fa - tj-i + a}"- 1 - 6 \ {tj _ 1<rj } 
+C [tj - tj-! + aY* l {r . <t ._ l} . 



3. Term A 7 _i7: 



— 1 _! 

fe-ij-i+a] 2 1 {t,- 1 <t J } - [tj -rj-i +a] 2 l {rj _ 1<tj} 
< C [tj-! - rj -i]*~ H+e [t, - tj-! + a] H_1 ~ £ l {Tj _ 1<t] _ 1<t]} . 
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4. Term Aj_iA 3 -: 

= - tj-i + <A 2 J-^-^t,} - fa - tj-i + a \ 2 1 {t 3 -i<r J } 

- - rj-i + af 2 l{r,-_i<t J } + fa - rj-i + af 2 l {r ._ 1<r . } 
< C[ tj - r^ H+£ - ri_i]'- ff4 * [r, - + a] 2ff -i" 2£ 1 



{tj_i<rj<t :j } 



_ 1 

+C [Tj - r --i + a] 2 l {rj _ 1<rj<tj _ 1<tj} . 



l {r J <t J ^ 1 <t j } 



If we replace the constant a by Sj — and we treat the the term Sj — Sj-i 
in the same way, using the inequality 

(a + by a < o _ *6 - *, 

we obtain the same estimates as if we had started with 

/ f / fe-ti-i] 4 n if H(*i.' , j ) dr, 

* V fc j =i / 

instead of (|4. 13|) . As a consequence, it suffices to control the following integral 

2 



L (lp flt ' t]dt ) *• 



(4.14) 



where a, b 6 {0, 1}, and A,- has one of the following forms 



4° 


= fa 


-tj 




1 {t J - 1 <r j }, 




4i 


= [tj 


~ r j 






L {ti-i<rj} 


40,1 


= [tj 


~tj 


1 

-l] 4 




1 {'"i<*i-i}' 


.1,0 


= [tj- 


-1 — 


r 3-l\ 




iH-f-e -, 
lj - L {r J _i<t ;j _i<t J }i 


4 1,1 


= [tj 


~ r 3 


-1 +E 




[ r j-tj-l\ l{t J -_ 1 <r 3 <t 3 } 


A 1 ' 1 


= fe- 


-1 — 


r 3-l\ 


— tj- 




A 1 ' 1 


= h 


- r 3 


-l]" 4 




fe-1 _r J-l] ff 2 l{r J _ 1 <r 3 <t 3 -_i<t J }, 



and with the convention that any term of the form Aj or Aa' must be followed 
by A ' or A -' 1 and any term of the form A/ or A 1 .' must be followed by 
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Aj' or Aj' . It is not difficult to check that the integral (|4.14[) is finite. For 
instance, for a product of the form A^'^A^'l we get 



[r,-_i - 1 



{rj-i<tj-i<rj<tj} 
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U 1-1+e r , i2ff-#-2e 



'■j-U 



x [tj - r ] 

h-i - tj-2] 



-l+e 



(7f 



i'-l 



l2.ff-4-e 



-l+e 



and the integral in the variable rj of the square of this expression will be finite 
because AH — | — 2e > — 1. 

So, we have proved that sup e E(j£) < oo for all fc. Notice that all these 
moments are positive. It holds that lim E ^ E(J £ J$) exists, and this implies the 
convergence in L 2 , and also in L p , for all p > 2. ■ 

On the other hand, if the initial condition of Equation (ll.ip is a constant 
K, then for all n > 1 we have 



(n\) 2 \\f n (-,t,x)f Hr =K 2 E 



{K* H ) Sq(B 1 — B 2 )(r, r)dr 



< oo, 



provided H £ (|, |). In fact, by Lemma FOl we have 



(n!) 2 ||/ n (-,i,x)||^ r = K 



'[0,t]3» 



t=i 

-2 



'[M 2 ™ i=l i=l 



and it suffices to apply the above proposition. 

However, we do not know the rate of convergence of the sequence ||/n(-,t, 
as n tends to infinity, and for this reason we are not able to show the existence 
of a solution to Equation (jl.ip in this case. 



5 Moments of the solution 

In this section we introduce an approximation of the Gaussian noise W H by 
means of an approximation of the identity. In the space variable we choose the 
heat kernel to define this approximation and in the time variable we choose a 
rectangular kernel. In this way, for any e > and 6 > we set 

W?;!=( f Mt-s)p s (x-y)dW» y , (5.1) 
Jo Jm d 
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where 

<Ps{t) = ^1[ ,5](0- 

Now we consider the approximation of Equation (jTTTJ) defined by 

We recall that the Wick product u t ' S x oW^ x is well defined as a square integrable 
random variable provided the random variable u e t ' x belongs to the space D 1,2 
(see (|2.9p l. and in this case we have 

<] S y o Wf* = J' J Ms - r) Pe (y - zK : l6W r H z . (5.3) 
The mild or evolution version of Equation Ij5.2|) will be 



u\% = p t u (y) + / pt- s ^-yK%oW^dsdy. (5.4) 

JO JR d 

Substituting (|5.3p into (|5.4p . and formally applying Fubini's theorem yields 
u ti = Ptu (y)+ / / / Pt- S (x - y)ip s (s - r) Pe (y - z)u E s ' S y dsdy ) SW^ Z . 

JO JR d \JQ JR d / 

(5.5) 

This leads to the following definition. 

Definition 5.1 An adapted random field u e,s — {u £ tx ,t >0,i£ R d } is a mild 
solution to Equation 115.2$ if for each (r,z) G R+ x M. d the integral 

Y r,f =11 Pt-s{x - y)tps{s - r)p e (y - z)u e s s y dsdy 
Jo Jr 

exists and Y t,x is a Skorohod integrable process such that Ii5. 5\) holds for each 
(t,x). 

The above definition is equivalent to saying that u £ t ' x G L 2 (Q), and for any 
random variable F e D 1,2 , we have 

E(Fuf x ) = E{F)p t u {y) 

I pt- s (x - y)tps(s - -)Pe(y - ^ul^dsdy) ,DF) . 

JR d J I Hi 



Our aim is to construct a solution of Equation (|5.2p using a suitable version 
of Feynman-Kac's formula. Suppose that B = {B t ,t > 0} is a d-dimensional 
Brownian motion starting at 0, independent of W . Set 

! W^ x+B ds = [ [ [ Mt-s-r)Pe(B s +x-y)dW r H y ds 

JO Jo Jo JR d 

1 f Ap s y dW r H y , 

JR d 
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where 

Ki= I Vs(t- s - r)p e (B s + x - y)ds. (5.6) 
Jo 



Define 



1 



u e t % = E* \ u (x + B t ) exp (jf A s r]y dW r " y - -a £,a I I , (5.7) 
where a^ 5 = |U £ ' 5 ||!, . 

II l\ ri d 

Proposition 5.2 T/ie random field u E 'f given by f<5. 9f ) is a solution to Equation 



Proof The proof is based on the notion of 5 transform from white noise 
analysis (see [5]). For any element tp £ Hi we dehne 



S, 



't, x (tp) = E (u E t ^F v J , 



where 

F v ,=exp (V H fo>) " 5 IMI* d 



From (|5.7p we have 



S t , x (<p) = E[u (x + B t )exp[W H (A*> s + ip)--a e ' s -~\\ipf nd 



E(ua(x + B t )exp((A^ s ,ip) Hd ^ 



= E uo 



)(x + B t )exp ^ (<p s (t - s - -)p E {B s +x - •), tp) Ud ds^j 



By the classical Feynman-Kac's formula, St, x {*p) satisfies the heat equation with 
potential V(t,x) — (ips(t — ■)p e (x — , that is, 

= ^S t , x (tp) + S t , x (cp) (lfis(t - -)p £ (x - •), ¥>} Hd ■ 

As a consequence, 



St,x(<p) =ptuo(x) + / / p t - s (x - y)S s , y (ip) (ips(s - -)p s {y - ■),ip) Hi dsdy. 



Notice that DF V = ipF v . Hence, for any exponential random variable of this 
form we have 

E (v*t£ F <p) = Pt u a{x) + / p t -s(x - y)E (u s t ]l (ip s {s - -)p s (y - -),DF v ) n ^ , 
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and we conclude by the duality relationship between the Skorohod integral and 
the derivative operator. ■ 

The next theorem says that the random variables u £ t ' S x have moments of all 
orders, uniformly bounded in e and S, and converge to the solution to Equation 
as 5 and e tend to zero. Moreover, it provides an expression for the 
moments of the solution to Equation 

Theorem 5.3 Suppose that H > ^ and d = 1. Then, for any integer k > 1 we 
have 



sup E 

e,S 







k~ 













< oo, 



(5.8) 



and the limit lim^o hm^o u \ 'x exists in IP , for all p > 1, and it coincides with 
the solution Ut }X of Equation Furthermore, ifU B (t,x) = Ylj=i u o(x+B{), 

where B 3 are independent d- dimensional Brownian motions, we have for any 
k>2 



E [<J = E 



B 



U*(t,x)exp I So{B\-Bi)d, 

\i<i Jo 



if H — \, and 



E [< J = E 1 



U*(t,x)e*j>{ f I 



4>(s,r)5 (Bi-Bi)dsdr 



(5.9) 



(5.10) 



ifH > |. 

In the case d — 2, for any integer k > 2 there exists to(k) > such that 
for all t < to(k) A5.8\) holds. If t < to(M) for some M > 3 then the limit 
lim e ^o lim,5|o u\' x exists in LP for all 2 < p < M , and it coincides with the 
solution Ut x 

of Equation [TJ]) . Moreover, l5A(f\) holds for all 1 < k < M - 1. 

Proof Fix an integer k > 2. Suppose that B l = {B\,t > 0}, i = 1, . . . , k are 
independent d-dimensional standard Brownian motions starting at 0, indepen- 
dent of W H . Then, using (|5.7p we have 



E 




uq(x + B 3 t ) exp 



Ai^dW! 1 



r 



where Apy'- 83 and ct e ' S ' B are computed using the Brownian motion B 3 . There- 
fore, 



E 











E B 


exp 










3 = 1 



,e,S.B J 



Y[u (x + B{' 



= E 



exp 



£ / A c,5,B' )Ae ,5,B>' 



i<3 



Y[u (x + B}) 

3 = 1 
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That is, the correction term ^a e ' S in (|5.7p due to the Wick product produces a 

cancellation of the diagonal elements in the square norm of Ylj=i A £ ' S ' B1 . The 

next step is to compute the scalar product (^A e - S ' B \ A e - S ' B3 ^j for i ^ j. We 

consider two cases. 

Case 1. Suppose hrst that H = \ and d = 1. In this case we have 

'A^ B \ A^ BJ ) = f [ [ f Vs (t-8 X - r) Pe (Bl i+ x-y) 
' ni j Jo Jo 

Xip s (t - s 2 - r)p e (B J S2 + x - y)dsids 2 drdy 



<ps(t - si - r)tps(t - s 2 - r) 

10 Jo Jo 
xp2e{B l Si - Bl 2 )d Sl ds 2 dr. 

We have 

/ (ps(t - si - r)(fs(t - s 2 - r)dr 
Jo 

= <T 2 [(< - si) A (i - s 2 ) - (t - Sl - S)+ V (t - s 2 - 5)+] + 
= Vs(si,s 2 ). 

It it easy to check that rjs is a a symmetric function on [0,t] 2 such that for any 
continuous function g on [0,£] 2 , 



lim / / rjs(si,s 2 )g(s 1 ,s 2 )ds 1 ds 2 
s l° Jo Jo 



g(s,s)ds. 



As a consequence the following limit holds almost surely 

\irn(A^ B \A^ B3 \ = f p 2e {^ - B{)ds, 

<HU \ / Hi Jo 

and by the properties of the local time of the one-dimensional Brownian motion 
we obtain that, almost surely. 



limlim(A e ^ B \A e ^ BJ 

elO SIO \ I Hx 



S (Bi-Bi)ds. 



The function rjs satisfies 

sup / r) S (s,r)ds < 1, 

0<r<t JO 

and, as a consequence, the estimate Q3.8P implies that for all A > 



supi? 

e,5 



Aexp (A 



\e,S,B' j^e,S,B 



Hi 



< OO. 
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Hence (|5.8p holds and lim e ^o hm^o u\ 'f := Wt i2; exists in L p , for all p > 1. 
Moreover, E{vf x ) equals to the right-hand side of Equation (|5.9p . Finally, 
Equation f|5.5[) and the duality relationship (|2.8[) imply that for any random 
variable F € D 1,2 with zero mean we have 



and letting <5 and e tend to zero we get 



Hi, 



E{F "" )=E {{ DF 'fJML 



Hi 



which implies that the process v is the solution of Equation (jl.ip . and by the 
uniqueness Vt <x = Ut, x - 

Case 2. Consider now the case H > h and d = 2. We have 

A £AB U MlBi ) = / fff [ t Mt-s 1 -r 1 )p e (Bl 1 +x-y) 

X(p s (t - s 2 - r 2 )p e (B J S2 + x-y)ds 1 ds 2 <j)(r 1 ,r 2 )dr l dr 2 dy 
rt ft ft ft 

tps(t - si - ri)ips(t -82- r 2 ) 

/o Jo Jo Jo 
xp 2s (Bl l - Bl 2 )ds 1 ds 2 (j)(ri 1 r 2 )dridr 2 . 

This scalar product can be written in the following form 

A e,6,B* }A eAB'\ = f f ^ - S 2 )p 2e {B\ i ~ BQds^, 

1 H d Jq Jq 

where 

r}s{si,s 2 )= / / ips(t-s 1 -r 1 )ip s (t-S2-r 2 ) (f){ri,r2)dridr2- (5.11) 
Jo Jo 

We claim that there exists a constant 7 such that 

m(8i,s 2 ) <l\si~s 2 \ 2H - 2 . (5.12) 
In fact, if \s 2 — s\\ = s we have 

s+8 



V8(si,s 2 ) < H(2H-1)S- 2 f f \u-v\ 2H - 2 

Js-8 Js 



dudv 



_L [ (s + 5 fH _ {s _ S) 2H _ 2g 2 H] 



Then 

HS- 2 [ S+S iy 2 "- 1 -(y- 6) 211 - 1 ) dy < H5 2H - 2 < H2 2 - 2H s 2B 2 
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if s < 25. On the other hand, if s > 26, we have 



< H(2H- l){s~5) 2H -' 2 

< H(2H-1)2 



2-2H g 2H-2 



It it easy to check that for any continuous function g on [0,i] 2 , 



/ / r] S (s 1 ,S2)g(si,s 2 )ds 1 ds2 = / / 4>(si, s 2 )g{si, s 2 )ds 1 ds 2 - 
lo Jo Jo Jo 

As a consequence the following limit holds almost surely 



lim 

.510 



lim\im(A e > S ' B \A £ > S ' B3 

ej.0 .510 \ / H d 



(sx,s 2 )&q(B\^ -B J S2 )ds 1 ds 2 . 



From (|5.12[) and the estimate (|3.13p we get 

sup£; B exp^A (a^ b \A s ^ b3 



< 00, 



(5.13) 



if A < Ao(t), where \o{t) is defined in (|3.12p with gammar replaced by 7. 

Hence, for any integer k > 2, if t < to(k), where fc ^~^ = Xo(t (k)), then 
holds because 



E 



( u ti) 



< \\uo\\ k E B 



\ 2 \ / H d 



Finally, if t < t (M) and M > 3, the limit hm £ | um 5|o u t'x := v t,x exists in L p , 
for all 2 < p < M and it is equal to the right-hand side of Equation (|5.10D . As 
in the case H — i we show that vt, x = u t , x - ■ 



6 Pathwise heat equation 

In this section we consider the one-dimensional stochastic partial differential 
equation 

^ = ±A U + «W£, (6.1) 

where the product between the solution u and the noise W B X is now an ordinary 
product. We first introduce a notion of solution using the Stratonovich integral 
and a weak formulation of the mild solution. Given a random field v = {v t x , t > 

0,a; e K} such that Jq J R \vt,x\ dxdt < 00 a.s. for all T > 0, the Stratonovich 
integral 

f f v tiX dW t H x 
Jo Jr 
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is defined as the following limit in probability if it exists 

i-T 

vt, x Wt'£dxdt, 



lim lim 

e|0 610 J 



where W^' x is the approximation of the noise W H introduced in (|5.1 



Definition 6.1 A random field u — {u tjX ,t >0,i£ M} is a weak solution to 
Equation \6.1\) if for any C°° function ip with compact support on R, we have 



ut, x <p(x)dx = / uo(x)(f(x)dx+ / / u s , x ( p"{x)dxds-\- / / u s ,xf(x)dW s x . 

Jm. JO J~R JO Js, 

Consider the approximating stochastic heat equation 

(6.2) 



Theorem 6.2 Suppose that H > 4. For any p > 2, the limit 

Umlimwf'f = u t x 

ej.0 510 ' 

exists in LP, and defines a weak solution to Equation in the sense of 

Definition \6.1\ Furthermore, for any positive integer k 



V / / (/>(s 1 ,s 2 )S(B i Sx -Bi 2 )ds 1 ds 2 
l0=l Jo Jo 



where U^fax) has been defined in Theorem \5.3\ ). 
Proof By Feynman-Kac's formula we can write 

<i = E B Su (x + B t ) exp (J J Ap s y dW r H y 

where Ap y has been defined in (|5.6p . We will first show that for all k > 1 



(6.3) 



supE 

6,e 







fe" 




e,6 











< oo. 



(6.4) 



Suppose that B % = {B\,t > 0},i = 1, . . . , k are independent standard Brownian 
motions starting at 0, independent of W H .Then, we have, as in the proof of 
Theorem [ 



E 



= E 



exp 



\e,6,B 1 j^e,6,B 



'■J = l 



Hi 



U?(t,x) 
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Notice that 



A e,SB* iA e,6B'\ = f f ^ ^^(^ _ B^ds^ 
1 Wi Jo Jo 

where r]s(si,S2) satisfies (j5.12j) . As a consequence, the inequalities (|3.8[) and 
(|3.15[) imply that for all A > 0, and all i,j we have 



sup£ expX(A^ 5B \A^ 5B3 

e,8 \ \ 'Hi 



< CO. 



Thus, holds, and 



lim lira E 

e{0 510 



E exp 



U B (t,x)exp [\ V / / - Bl)d Sl ds 2 

\ 1 i, i= i -/o 



In a similar way we can show that the limit lim E) e'j.o lirri<5,<5'lo E (u £ t \ 



8 e',6' 



ists. Therefore, the iterated limit lim e m limaio E 
Finally we need to show that 



e,8 



exists in L 2 . 



lim lim 

ej.0 510 



f f u SiX cp(x)dW s H x - f f //.;•;,;■:,• ill . ",</>,/,•) = 0, 

o Jr Jo Js. / 



in probability. We know that f Q L ul' x (p(x)W^ x dsdx converges in L 2 to some 
random variable G. Hence, if 



B 



e,8 



u E s ' x — u SyX ) ip(x)W sx dsdx 



(6.5) 



converges in L 2 to zero, u SiX cp(x) will be Stratonovich integrable and 



u s ,x l p(x)W s x dsdx = G. 



The convergence to zero of (|6.5p is done as follows. First we remark that B Ei s 
6(<j) £ ' S ), where 



(u e s ' S x — u SjX ) tp(x)(fis{s — r)Pe( x ~ z)dsdx. 



Then, from the properties of the divergence operator, it suffices to show that 



lim lim \\D(jf' S 



-:\0 510 



(6.6) 
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ke,6\ 



= 0. On the other hand, 



It is clear that lim^o lim^o E \ \\ Y \\ Hi 

M / ' (D (u e s f x ) - D (u s<x )) <p(x)ip s (s - r)p s (x - z)dsdx, 



D 



and 



D(u 



E 



uq (x + B t ) exp 



Then, as before we can show that 



lim lim E\ (D 

e,£'|0(5.<5'J,0 



= E 



(x + B})u (x + B 2 t ) exp ( I I ^s 2 )5 Q {B l Si - Bi 2 )d Sl ds 2 
\ iJ=1 Jo Jo 



[ [ ^( Sl) S2Mo« -Bl)d Sl ds 2 
Jo Jo 



This implies that u e s ' x converges in the space ID 1 ' 2 to u SlX as 5 J, and e J. 0. 
Actually, the limit is in the norm of the space ID) 1 ' 2 (Hi). Then, (|6.6p follows 
easily. ■ 

Since the solution is square integrable it admits a Wiener-Ito chaos expan- 
sion. The explicit form of the Wiener chaos coefficients are given below. 



Theorem 6.3 The solution to i6.1\) is given by 

oo 

"4,0! = / ]ln(fn(;t,x)) 



(6.7) 



n=0 



where 



fn it\ : X\ , . . . : t n , X n , t, x) 



= E L 



u (x + B t ) exp ^~ J J 4>(s 1 ,s 2 )5 (B Sl - B S2 )ds 1 ds 2 
x8 (B tl +x - £1) • ■•5o(-B tn + x - x n )\ . (6.8) 
Proof From the Feynman-Kac formula it follows that 

E B (u (x + B t ) exp (J A^dW^Vj 

E B \u {x + B t ) exp exp Qf* J^dW r H y - \\\A^ s f ni 

OC 

^Inif^^x)), 
n=0 



E,S 
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where 



u (x + B*) exp A^ X1 • • • ^ 



Letting S and e go to 0, we obtain the chaos expansion of ut, x - ■ 

Consider the stochastic partial differential equation (|6.ip and its approxima- 
tion (|6.2[) . The initial condition is uo(x). We shall study the strict positivity of 
the solution. In particular we shall show that E \\ut{x)\~ p ] < oo. 

Theorem 6.4 Let H > 3/4. If E (\u (B t )\) > 0, then for any < p < oo, we 

<oo (6.9) 

and moreover, 



E [\u t (x)\-P] < (E\u (x + Bt)])-"- 1 E B 



u (x + B t )\ 
exp J J S(B Sl - B S2 )4>(si,s 2 )dsids 2 



(6.10) 



Proof Denote k p = (E B (\uq(x + B t )\)) P . Then, Jensen's inequality ap- 
plied to the equality u £ t ,S x = E B ^uq{x + B t ) exp ([* J R ApydWjFy) j implies 
that t 

\uti\~P < Kp E B [ \u (x + B t )\ exp (-pj o jf ^dW* 



Therefore 



e,<5i 



< k p S a ^ |u (a; + B t )|E 



= K p E 



\u (x + B t )\E 



i'X|) j — p 

2 



Ap 5 y dW B y 



and we can conclude as in the proof of Theorem [ 

Using the theory of rough path analysis (see [8]) and p- variation estimates, 
Gubinelli, Lejay and Tindel [4] have proved that for H > |, the equation 

§f -1*. 

had a unique mild solution up to a random explosion time T > 0, provided 
a e Cj(R). In this sense, the restriction H > |, that we found in the case 
u(a:) = a; is natural, and in this particular case, using chaos expansion and 
Fcynman-Kac's formula we have been able to show the existence of a solution 
for all times. 
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